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We propose a generalization of multifractal analysis that is applicable to the critical regime of 
the Anderson localization-delocalization transition. The approach reveals that the behavior of the 
probability distribution of wavefunction amplitudes is sufficient to characterize the transition. In 
combination with finite-size scaling, this formalism permits the critical parameters to be estimated 
without the need for conductance or other transport measurements. Applying this method to high- 
precision data for wavefunction statistics obtained by exact diagonalization of the three-dimensional 
Anderson model, we estimate the critical exponent v = 1.58 ± 0.03. 
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The statistical analysis of spatial probability and den- 
sity fluctuations, which has a distinguished history [l|, 
has recently received new impetus. Scanning-tunnelling 
spectroscopy now allows the direct measurement of the 
spatial variation of charge densities y. Dramatic ad- 
vances in cold atom physics are stimulating the study of 
Anderson localization in Bose-Einstein condensates by 
imaging of the atomic densities. This permits the obser- 
vation of the exponential decay of the wavefunctions and 
direct measurement of localization lengths [3| . Anderson- 
type transitions can now be investigated experimentally 
in quasi-periodic disorder potentials [J] and in cold-atom 
realizations of the kicked-rotor [5| . Similarly, the spatial 
localization of light [6| has recently been studied in nano 
devices with slow- wave structures [7|. The fundamental 
tool to characterize density fluctuations and the scale in- 
variance of spatial distributions at critical points is mul- 
tifractal analysis (MFA). Recently, a predicted symmetry 
of the multifractal spectrum at the Anderson transition 
[8| has been confirmed by experimental studies of vibra- 
tions in elastic networks [9|. MFA has also furnished 
insights into the theoretical foundations of the quantum 
Hall transition [lOJ. However, the application of MFA is 
restricted to the critical point, where the relevant prob- 
ability distributions are truly multifractal [ll|, |l2|. Up 
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to now, this has meant that additional computer simu- 
lations or experiments must be performed in advance to 
locate the critical point precisely, since any error here 
will adversely affect the results. 

In this Letter we propose a generalized MFA that is ap- 
plicable throughout the critical regime and not just at the 
critical point. Our approach is motivated by the behavior 
of the probability density function (PDF) of the wave- 
function intensities, ■p(a; W, L, A). Here, a = In/ifc/lnA 

with /j,fe = J2i=i IV'iP the summed wavefunction proba- 
bility in the k-th cubic box of linear size ^ in a lattice 
with volume L^, and A = i/L. We have applied our 
method to the three-dimensional Anderson model with 




FIG. 1: (color online) Evolution of the wavefunction ampli- 
tude distributions P{a; W, L) as a function of disorder W 
across the Anderson transition, at fixed A = 0.1 for two system 
sizes L. Each distribution has been computed with 10^ wave- 
functions. The data points (•) and solid lines on the bottom 
plane mark the trajectories of the maximum 5m. For clarity, 
distributions are shown a,t W = 15, 16.6 and 18.0 only. 



box-distributed site energies of width W. We have cal- 
culated more than 1.5 million uncorrelated wavefunctions 
by exact diagonalization of system sizes up to 100'^ [l3|. 
We find that the parameter dependence of P{a; W, L, A), 
as displayed in Fig. [H is sufficient to characterize the 
Anderson transition. For fixed A, the distribution be- 
comes scale invariant at the critical point and away from 
the transition its maximum, Sm, exhibits finite size scal- 
ing (FSS) behavior: 5m shifts in opposite directions in 
the different phases at a rate which depends on L. This 
provides an alternative way to estimate the critical pa- 
rameters of the transition that is not based on transport 
properties such as the conductance. We believe that our 
approach is particularly valuable in experiments where 
the PDF of wavefunction amplitudes is accessible, e.g.. 



A V 



Wc 



Nd Np X P nonimgnir, 



T2 0.1 1.58(52,66) 16.57(50,61) 153 13 151 0.2 5 13 

T2 0.2 1.59(57,61) 16.56(52,58) 153 12 158 0.2 6 2 

T3 0.1 1.62(57,66) 16.56(42,61) 153 10 133 0.7 5 10 

5m 0.1 1.56(54,59) 16.53(49,55) 153 10 131 0.7 3 2 10 

5i 0.1 1.61(58,64) 16.56(52,59) 81 7 89 0.1 2 10 

TABLE I: The estimates of the critical parameters v and Wc, 
together with 95% confidence intervals, from one-parameter 
FSS. The number of data is Nd, the number of parameters 
is Np, x^ is the value of the chi-squared statistic for the best 
fit, and p is the goodness of fit probability. The orders of the 
expansions are specified in the last column. The system sizes 
are L G [20, 100], and the range of disorder is M^ € [15, 18] 
except for 5i where W £ [16, 17]. The data uncertainty is on 
average 0.6% for rs, 0.3% for T2 and 0.1% for 5i and 5m. 
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FIG. 2: (color online) Two-parameter finite-size scaling result 
for 5m (A) at energy E = when fitted according to Eq. ([2]) 
indicated by the shaded surface. Data (•, o) and fit parame- 
ters are as described in the second row and caption of Table 
im The black lines on the surface highlight the different values 
of A [symbol (o) highfights A = 0.1]. The estimated ao = 4.09 
corresponds to the extrapolation of 5m (A) as A — >■ 0. 



through LDOS measurements using STM techniques jj] 
or in ultracold Bose/Fermi gases in disordered optical lat- 
tices [sl, IJ] . For the Anderson model our estimates of Wc 
and v in Table H] are in excellent agreement with previ- 
ous transfer matrix results [14;], which resolves the long- 
standing issue of systematically smaller exponents found 
in previous diagonalization studies 15] . In addition, the 
size dependence of the PDF on L and X (or, equivalently, 
i) suggests the possibility of two-parameter FSS. This is 
illustrated in Fig. [5] where scaling is shown as a function 
of both L/^ and A, with ^ the correlation/localisation 
length. 

Our generalization of MFA starts by considering the 



g-moments of the wavefunctions defined as Rq = J^k Mfej 
where the sum runs over the (L/i)'^ boxes of linear size 
£. At criticality and in the thermodynamic limit (L — ^ 
cxj), due to the multifractal nature of the states, the only 
relevant parameter is A = i/L and the moments scale 
as (Rq) ^ X^'' . Here, the brackets denote an average 

over disorder. Away from the transition, however, the 
moments depend on £, L and the disorder W . It follows 
from scaling arguments that close to the transition Rq — 
Rq{L/i, A) = RqiL/^,e/0 [H, with C = aW) diverging 
at the critical point Wc as ^ oc \W — Wc\ 
criticality we can write {Rq){W,L,X) = X'^ 
which can be rearranged as follows, 
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TZqiL/^X), 



TqiW,L,X) 



g(g-l) 
In A 



rq{L/^,x). 



(1) 



Here, Tq is related to the original TZq and we have defined 
a generalized mass exponent as Tq{W, L, X) = \n{Rq)/ In X 
which becomes the usual Tq at Wc and in the limit A ^■ 
0. The factor q{q — 1) has been explicitly included to 
satisfy tq = tq = —d and ri = ti = 0. From Eq. ([1]) 
it is straightforward to obtain the scaling law for the 
singularity strengths 3g = d/Fq/dq, 

aq{W,L,X)^aq + :^Aq{L/^,X), (2) 

where the second term on the rhs will be non-zero for 
all q values, and the generalized exponents are defined as 
aq{W,L,X) = (X]fcMfelnAife>/((-R<?)ln^)- Consequently, 
we can define a W, L and A dependent generalized sin- 
gularity spectrum fq = qctq — Tq , obeying 



q 



f^{W,L,X) = fq + :^TqiL/tX) 



(3) 



Eqs. (IT])-([3]) suggest a wide range of generalized expo- 
nents that can be used to perform FSS and obtain Wc 
and v. In addition, the scale invariant multifractal expo- 
nents Tq , aq and fq at the critical point can be estimated 
from the same FSS study without the need to know Wc 
beforehand. Moreover, the use of different moments of 
the wavefunctions provides a test of the stability of the 
estimates for the critical parameters, as these should be 
g-independent. 

The generalized multifractal spectrum ([3]) is re- 
lated to the PDF of the wavefunction amplitudes as 
V{a; W, L, A) ex xd-nS;W,L,x)_ ^g ^g approach the ther- 
modynamic limit (A — > 0) at the critical point, this be- 
comes the usual relation Vx{a) ex A''"-'"'-"' |l7|. As shown 
in Fig.[Tl for fixed A, the PDF becomes scale invariant at 
the transition. The generalized exponents can also be cal- 
culated from the distribution Via'jW, L, X), which may 
be useful when the wavefunctions cannot be probed indi- 
vidually and only partial information about the PDF is 
accessible. For g = 0, we have 5o = {^'^ J2k ^^ A^fc/ In A) ~ 
{a), which corresponds to the mean value of the PDF. 



a) converges towards the position of the 

a) 



23| . While a^ and ao 



When L ^ oo 
maximum of the PDF a 
may differ quantitatively at finite L, they obey the same 
scaling law with the same critical parameters. Therefore 
the scaling of either the mean value or the position of the 
maximum of the PDF as a function of W^, L and £ may 
be used to estimate the critical parameters. We note that 
the scaling law ([2]) that we give here for our generalized 
multifractal exponents Sq and 5ni is different from the 
scaling laws suggested in the past 
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We first present results for standard one-parameter 
(L/0 FSS at fixed A values [13, El- Let Vx[W,L) de- 
note either rg(W,L, A) or aq{W,L,X). We introduce a 
set of fit functions which include two kinds of correc- 
tions to scaling, (i) nonlinearities of the W dependence 
of the scaling variables and (ii) an irrelevant scaling vari- 
able that accounts for a shift of the disorder value at 
which the T\[W,L) curves cross. We use T\{W^L) = 
Q{qL^/^ .rjL^^y^), where Q denotes the rhs of either Eq. 
(P) or ^ and g and 77 are the relevant and irrelevant 
scaling fields, respectively. The function Tx{W,L) is ex- 
panded to first order in the irrelevant scaling variable 
as Vx{W,L) = GoigL^'") + 'qL-\y\gi{gL^/''), and subse- 
quently Qs = Y^^'Lo'^skQ^L^^" ■ The fields g and -q are 
expanded in terms of w = iWc — W)/Wc up to order nig 
and m^, respectively, such that g{w) = X]m=i ^m'"'™' 
Tl{w) = J2m=a'^rnW™i with bi — cq = 1. The expansions 
of the fit functions are truncated at orders no, ni, nig, r/i,,. 
The orders of these expansions should be kept as low as 
possible, while giving an acceptable goodness of fit proba- 
blity p. We emphasize that the results of the FSS analysis 
are valid only if the goodness-of-fit is acceptable. 

In Fig. [3] (top) we show a fit for T2 data at A = 0.1. As 
L increases, the generalised exponent T2 approaches the 
metallic and insulating limits correponding to d{q—l) for 
W <Wc and for H^ > Wc, respectively. We emphasize 
that the value of T2 at the critical point is A-dependent, 
and it is only in the limit A ^^ that T2(A) converges 
towards the scale invariant T2. A similar behavior is ob- 
served in Fig. [3] (bottom) for ctm, the position of the 
maximum of the PDF. In this case ImiL^ao 5m = d in 
the metallic side and limL^cx) Sm = 00 in the insulating 
phase. The trajectories of Sm as a function of disorder 
for different L are also shown on the bottom plane of 
Fig. [H The PDF V{a; W, L) was obtained from the nu- 
merical histogram of wavefunction intensities [l7[. The 
position of its maximum was estimated by fitting V{a) 
to A'''^"^ where h{a) is a polynomial, which allows for the 
non-symmetric and non-Gaussian nature of the distribu- 
tion jl7|] . The precision of the 5m data was determined 
by performing the fit 100 times on independent distribu- 
tions obtained from subsets of 100 states each for every 
set iW,L). 

In Table m we show representative results for Wc and i^ 
from 3 fits of Tq and 2 fits for 5q at various A and q values. 
The analysis is based on a total of 1, 530, 000 wavefunc- 




FIG. 3: Plot of T2 (top) and Sm (bottom) at E — and 
A = 0.1 as a function of disorder at various system sizes L £ 
[20, 100]. The errorbars denote standard deviations obtained 
from averaging over disorder. The lines are plotted according 
to Eq. ([T]) and ^, respectively, for the fit parameters as in the 
first and fourth rows of Tableland with irrelevant exponents 
y = -1.8±0.4 and y = -1.65±0.12. The vertical dashed lines 
show the estimated Wc and its confidence interval is indicated 
by the grey region. The insets show the FSS functions in an 
enlarged view of the critical region. 



tions generated at energy E = for system sizes L^ in 
the range 20'^ to 100^, and disorder values W from 15 to 
18, where for each pair {W, L) we average over 10* inde- 
pendent states. Table HI shows that all fits give estimates 
of v which (i) are consistent with each other, (ii) agree 
with the transfer-matrix-method results v = 1.57(55, 59) 
14J and v — 1.62(55,69) [l8| and (in) are significantly 
larger than 1 and, within the accuracy, different from 
1.5 ^]. Our results also agree with the estimation of v 
obtained from the quantum kicked rotor which was re- 
cently realized experimentally using cold atoms [5'|. The 
large irrelevant shift of Wc seen in Fig.[3]is comparable to 
those observed for higher Lyapunov exponents [2l| . This 
might explain the variation in the estimated value of v 
from previous works based on exact diagonalization [15| , 
as only the use of very large system sizes can resolve this 
shift unambiguously. 



The two-parameter scaling suggested in Eqs. H]) and 
^ and shown in Fig. [5] is based on a scaling function of 



w. 



Nn N, 



P X 



p expansion 



T2 1.56(52,60) 16.57(55,59) 544 20 558 0.15 320130 

5m 1.56(55,58) 16.55(54,56) 544 16 494 0.85 220 210 

?i.i 1.60(55,64) 16.57(55,59) 224 11 232 0.17 210010 

5i 1.60(55,64) 16.56(54,58) 224 13 230 0.18 310010 

TABLE II: The estimates of Wc and v from two-parameter 
FSS. Labels are described in Table H] The last colum 
shows the orders of the expansion: nj,, n°, n^, nj, nig 
and m,,. The scale invariant multifractal exponents ob- 
tained from the fits are T2 = 1.21(20, 22), qq = 4.09(08, 10), 
n.i = 0.184(183,185) and ai = 1.93(92,94) respectively. 
The data used correspond to L G [20,100], i ^ 2 satisfy- 
ing 0.02 s; A s; 1/7 and W € [15, 18], except for the last two 
where W G [16.2,16.8]. 



the type 

^ ^ (4) 

for aq (similarly for Tq) where the dominant irrelevant 
scaling is determined by i [2j] . The functions Aq are ex- 
panded in their arguments and relevant /irrelevant fields, 
and the expansion is characterized by the indices n^ , n^ , 
n\, n], rUg and m^. The two-parameter scaling provides 
a simultaneous estimation of the critical parameters Wc: 
V and the scale invariant multifractal exponents r,, a^. 
As shown in Table HIl the estimated values for Wc and v 
from two-parameter FSS, using a large number of data 
for integer and non-integer g, are in agreement with those 
obtained at fixed A (Table IJ). The data involved in this 
analysis may differ in i but share the same L, hence there 
is a certain degree of correlation that may affect the fit 
[22| . However, the values of v in Table HIl are within the 
accuracy the same as those of the uncorrelated FSS in 
Table H] [23. In Fig. [2] we show the scaling for ai„. The 
shaded surface denotes Eq. ^ using the two-parameter 
scaling function ^ with the irrelevant correction sub- 
tracted, displayed as a function of L/^ and A, where the 
correlation length is given by ^ = \g\~^ ■ The scaling 
function exhibits an upper and a lower sheet populated 
by values of Sm corresponding to extended (W < Wc) 
and localised (W > Wc) states respectively. The merg- 
ing of the two sheets as ^ — >■ 00 at constant A determines 
the estimation of v. The additional extrapolation of the 
merging point as A -t- gives the scale invariant a^. 

In conclusion, we have proposed a generalisation of 
multifractal concepts such as mass exponents, singularity 
strengths and the multifractal spectrum that is applica- 
ble to the critical regime of the Anderson transition. The 
combination of the generalized MFA with FSS provides 
the critical parameters of the transition and enables MFA 
to be applied without knowing the exact position of the 
critical point in advance. We have tested our method on 
the Anderson model of an electron in a disordered system, 
and we estimate the critical exponent that describes the 



divergence of the localization length to be i^ = 1.58±0.03, 
in agreement with previous transfer matrix calculations 
[IJj . The method is applicable to other models with crit- 
ical fluctuations and to the wealth of such experimental 
data that is now becoming available [2-y, LZ| ■ 
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